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Captain  Kevin  Malloy 
AFOSR-NE 
Bldg.  210  Rm.  219 
Bolling  AFB 

Washington  D.C.,  20332-6448 


Dear  Kevin, 


Please  consider  this  my  progress  report  for  one  year  since  Nov.  1986.  I  enclose 
copies  of  two  publications. 

Very  briefly  our  effons  this  past  year  have  been  directed  toward  the  development  of 
an  ab-initio  formalism  to  study  grain  boundaries  in  solids. 

The  microscopic  characterization  of  grain  boundaries  in  polycrystalline 
semiconductors  has  become  essential  to  our  understanding  of  the  physical  properties  of 
these  technologically  important  materials.  From  a  theoretical  standpoint,  this  requires  a 
formal  quantum  mechanical  treatment  of  the  total  energy  of  the  grain  boundary  core  region 
which  should  properly  take  into  account  the  redistribution  of  electron  charge  density  and 
the  unrestricted  relaxation  of  the  ions.  First  principles  approaches  to  calculating  the  total 
energy  of  solids  have  begun  to  be  applied  to  grain  boundaries  in  semiconductors  and 
metals.  However,  the  methods  employed  in  these  studies  have  not  been  entirely  ab  initio  in 
the  sense  that  assumptions  regarding  the  starting  atomic  configuration  were  made  or  that  the 
potential  functions  were  semi-empirical..  Recently,  how-ever,  we  succeeded  in  rierforming  a 
completely  ab  initio  investigation^of  the  structure  of  a  grain  boundary  in  geranium.  This 
was  accomplished  using  a  quantum  molecular  dynamics  simulated  annealing  method.  This 
method  allows  global  minimisation  of  the  boundary  energy  to  be  achieved  w  ith  respect  to 
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all  electronic  and  ionic  structural  degrees  of  freedom  using  ab  initio  local  pseudopotentials. 
The  method  has  significant  advantages  in  computational  speed  and  storage  requirements 
over  traditional  total  energy  techniques,  especially  when  systems  of  low  symmetry  are 
involved  or  in  which  large  relaxations  take  place. 

Various  interesting  results  emerge  including  identification  of  two  low  energy  S=5 
(001)  twist  boundary  structures,  calculation  of  boundary  formation  energies,  determination 
of  boundary  volume  changes,  and  identification  of  novel  intrinsic  grain  boundary  defects 
such  as  four-fold  rings  of  bonds  and  overcoordinated  atoms. 

At  present,  we  are  performing  calculations  of  the  valence  electron  charge  density  in 
the  neighborhood  of  the  boundary,  in  order  to  determine  the  exact  nature  of  local 
interatomic  bonding.  We  are  also  engaged  in  mapping  out  the  total  energy  in  a  dense  grid 
of  translation  states  in  order  to  obtain  a  global  view  of  the  energy  surface. 

Sincerely, 

J.  D.  Joannopoulos 

Professor  of  Physics 
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Ab  Initio  Determination  of  the  Structure  of  a  Grain  Boundary  by  Simulated  Quenching 


M.  C.  Payne.''*’  P.  D.  Bristowe,  and  J.  D.  Joannopoulos 
Stassachusetls  Institute  of  Technology.  Canthridge.  Massachusetts  021 19 
(Received  5  January  1987) 


Results  of  the  rtrsl  completely  ab  initio  investigation  of  the  microscopic  structure  of  a  grain  boundary 
in  a  semiconductor  are  presented  By  use  of  the  moiccular-dynamics-simulated  annealing  method  for 
performing  total-energy  calculations  within  the  local-density- functional  and  pseudopotential  approxima¬ 
tions,  the  I  ”5  (001)  twist  boundary  is  germanium  is  studied.  .\  number  of  rotation-and-translation 
states  are  investigated  leading  to  a  prediction  for  the  structure  of  this  geometry.  Evidence  for  the  possi¬ 
ble  presence  of  novel  delects  and  glasslike  tunneling -mode  stales  at  grain  boundaries  is  presented. 


PXC'S  numbers  61  '0  Ng,  68  35  Ja,  68  35  .Md 

Internal  interfaces  in  a  solid,  such  as  grain  boundaries, 
often  control  the  mechanical,  chemical,  and  electrical 
properties  of  the  material.  It  is  essential,  therefore,  to 
have  an  understanding  of  their  structure  at  an  atomic 
level  Unfortunately,  the  structural  complexity  of  grain 
boundaries  makes  their  study  a  formidable  task  both  ex¬ 
perimentally  and  theoretically.  Recently,  experimental¬ 
ly,  new  electron-diffraction  and  imaging  techniques  have 
begun  to  yield  detailed  structural  information  about  a 
number  of  gram  boundaries  in  special  tilt  and  twist 
orientations.'*-  On  the  theoretical  side,  first-principles 
approaches  to  calculating  the  total  energy  of  solids  have 
begun  to  be  applied  to  interfaces. These  studies  have 
provided  important  information  on  the  distribution  of 
charge  density  at  the  boundary,  the  local  density  of 
states,  the  atomic  structure,  and  the  boundary  energy. 
However,  the  methods  employed  so  far  have  not  been  en¬ 
tirely  ah  imtio  in  the  sense  that  a  starting  atomic- 
configuration  has  been  assumed  based  on  geometric 
models  or  experimental  observations. 

In  this  paper  we  present  a  completely  ah  initio  investi¬ 
gation  of  the  structure  of  a  twist  grain  boundary  in  ger¬ 
manium  where  nothing  is  known  about  its  microscopic- 
structure  -X  method  is  employed  which  provides  for  the 
self-consistent  and  simultaneous  relaxation  of  the  atomic 
,ind  electronic  coordinates.  The  only  input  data  are  the 
atomic  number  and  mass  of  germanium. 

Specifically,  we  consider  the  I  “7  (001)  twist  bound¬ 
ary  m  the  diamond  structure,  where  I  is  the  inverse  den¬ 
sity  of  coincidence  sites.  Two  distinct  I  “5  boundary 
structures  with  the  same  coincidence  site  lattice  (CSIJ 
can  he  formed  hy  rotating  one  diamond  half-crystal  with 
respect  to  another  about  the  [00 1 1  direction  by  3f>  9"  or 
iiv  complementary  angle.  —53. 1”  "  Both  structural 
variants  arc  studied  in  the  present  work  and  are  desig¬ 
nated  15  and  15*  .-^s  an  example,  in  Fig.  1(a)  we  illus¬ 
trate  the  .itomic  positions  of  two  layers  on  either  side  of 
,in  ideal,  unrelaxed  15  boundary  In  this  case,  the  axis 
ot  rot.ition  passes  through  the  atom  at  the  center  of  the 
indicated  unit  cell  Lxperimentally,  it  is  not  possible  to 
constrain  the  position  of  this  axis  to  any  specific  transla¬ 
tion  si.ite  ,ind  .1  given  rotation  state  will  adopt  the 

I  -mH 


translation  state  of  the  lowest  energy.  Since  the  position 
of  ihis  rotation  axis  is  a  priori  unknown  to  within  a 
‘‘displacemenl-shift-complete''  vector'*  of  310  it 
becomes  a  challenge  to  theory  to  predict  the  appropriate 
lowest-energy  translation  state. 

It  is  not  possible,  however,  to  calculate  the  energies  of 
all  the  translation  states  within  a  displacement-shift- 
complete  lattice  unit  cell.  In  this  work  we  study  the  CSL 
configurations  for  both  rotation  states,  the  io  '-'O' 
translation  state  of  15  and  the  210)  translation  stale 
of  £5*.  These  particular  states  were  chosen  because  the 
structures  have  higher  symmetry  than  boundaries  pro¬ 
duced  by  arbitrary  translation  states  .Moreover,  in  the 
pair-potential  simulation'^  of  fee  twist  boundaries  in 
metals,  such  slates  appear  to  have  the  lowest  energies. 
Whether  this  is  true  for  semiconductors  is  certainly  not 
clear.  Nevertheless,  for  simplicity,  we  restrict  ourselves 
to  these  four  roiaiion-and-lranslation  states 

The  grain-boundary  structures  were  modeled  within 
the  supercell  approximation  which  imposes  periodicity 
normal  to  the  interface  bound.try  Experimentally,  it  is 
found  that  grain-boundary  '  widths”  in  semiconductors 
arc  extremely  narrow  (  :l  10  so  that  a  relatively 

small  number  of  atomic  layers  on  either  side  of  the  inter¬ 
face  can  accurately  model  the  system.  Our  supercell 
contains  twelve  layers  of  atoms  representing  two  six- 
layer  crystal  slabs  rotated  relative  to  each  other  by  either 
3b. O'*  or  —53.1°  .\U  the  atoms  in  the  supercell  were  al¬ 

lowed  to  relax  except  for  the  atoms  in  the  central  two 
layers  of  each  crystal  slab  which  were  confined  to  their 
crystallographic  positions  In  addition,  the  distance  be¬ 
tween  these  fixed  layers  was  allowed  to  vary  as  the  local 
gram-boundary  volume  changed  The  atomic  positions 
and  local  grain-boundary  volume  were  assumed  to  have 
reached  their  equilibrium  values  when  the  forces  had  be¬ 
come  smaller  than  0.1  eV/A. 

The  calculations  were  performed  with  use  of  the 
moiecular-dvnamics  simulated  quenching  method  for 
total-energy  calculations.'*''  The  exchange-correlation 
potential  was  calculated  within  local-density-functional 
theory  with  use  of  the  Perdew  and  Zunger  parametriza- 
tion.  The  electron-ton  interaction  was  determined  with 
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F'.G.  !  Atomic  positions  in  planes  normal  to  the  tOOt)  direction  for  two  layers  above  (open  circles)  and  below  (filled  circles)  a 
iwisi-grain  boundary  in  germanium,  (a)  Unrelaxed  15  CSL;  a  typical  unit  cell  is  indicated,  (b)  Relaxed  15  CSL.  (c)  Unrelaxed 
15  57 '210' state  (d)  Relaxed  15  ji  CIO’ state  (e)  Relaxed  15*  CSL  (f)  Relaxed  15*  15  (210' state 


use  of  the  Starkloff-Joannopoulos  ab  initio  local  pseudo¬ 
potential  approach.  The  molecular-dynamics  scheme 
is  particularly  suited  to  the  relaxing  of  systems  with 
many  degrees  of  freedom.  Briefly,  conventional  matrix 
diagonaltzation  is  avoided  and  replaced  by  simple  in¬ 
tegration  of  classical  equations  of  motion  for  the  elec¬ 
trons,  ions,  and  unit-cell  volume.  This  method  allows  the 
use  of  a  large  number  of  basis  functions  and  our  calcula¬ 
tions  were  performed  with  an  energy  cutoff  of  10  Ry, 
which  corresponds  to  about  4300  plane  waves.  Two 
Monkhorst-Pack^®  k  points  were  used  for  Brillouin-zone 
averaging.  No  symmetry  was  imposed  on  the  system  and 
the  symmetry-breaking  element  in  the  molecular- 
dynamics  method  allowed  the  systems  to  explore  regions 
of  phase  space  with  lower  symmetry  than  the  initial 
configurations. 

The  results  of  all  our  calculations  arc  presented  in 
Table  I  and  Fig.  1.  The  table  lists  the  total  energies  for 
the  four  twist  boundaries  considered  and  the  change  in 
grain-boundary  volume  after  relaxation.  The  geometric 
structures  of  the  rela.xed  boundaries  are  shown  in  Figs 
Kb),  1(d).  1(e),  and  1(f).  In  the  figures,  we  have  arbi¬ 
trarily  chosen  to  draw  bonds  between  all  pairs  of  atoms 


that  are  separated  by  distances  less  than  1.15  times  the 
bond  length  in  the  bulk.  In  the  extremely  distorted  envi¬ 
ronment  of  the  grain  boundary,  it  is  possible  that  atoms 
separated  by  this  distance  could  still  be  bonded  although 
a  detailed  investigation  of  the  electronic  charge  density  is 
required  to  test  this.  We  shall  present  a  detailed  study  of 
the  bonding  in  the  grain  boundaries  in  a  forthcoming 
publication. 

Figure  Kb)  shows  the  structure  obtained  after  allow¬ 
ing  the  (15  CSL)  grain  boundary  in  Fig.  Ka)  to  relax. 
The  relaxation  basicallv  involves  the  formation  of  a  di- 


TABLE  I  Gram-boundary  energies  (£)  and  local  volume 
changes  (Atl ). 


Energy  (£) 


Rotation 

36.9”  (15) 
36  9”  (15) 
-53.1  (15*) 
-53.1  (15*) 


CSL 

S'-210> 

CSL 

fo  '210  ’ 


(mJ  m  •) 

(eV) 

An  (A) 

1288 

5.87 

4-0  07 

851 

3  88 

-0  05 

1415 

0  45 

+0  07 

937 

4.27 

-0 ;; 

-  ^  ^  *  Ilk. 
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mer  bond  between  a  pair  of  atoms  in  the  layer  above  the 
interface.  We  note  that  it  appears  to  be  energetically 
unfavorable  for  new  bonds  to  form  on  both  sides  of  the 
interface.  The  distortion  leaves  six  atoms  per  unit  cell 
with  only  three  near  neighbors  and  the  structure  contains 
one  four-membcred  ring  per  unit  cell.  Fourfold  rings  of 
bonds  are  very  unusual  in  germanium  or  silicon  struc¬ 
tures  and  to  our  knowledge  have  never  been  observed  be¬ 
fore.  They  are,  however,  an  intrinsic  structural  element 
of  this  type  of  twist  boundary. 

In  Figs.  1(c)  and  1(d)  we  show  respectively  the  unre¬ 
laxed  and  relaxed  geometries  of  the  15  ^(210)  grain 
boundary.  The  relaxation  in  this  boundary  is  much 
larger  than  in  the  previous  case  and  some  atoms  move 
distances  greater  than  40%  of  the  bond  length  during  the 
relaxation.  .All  the  atoms  have  four  near  neighbors  after 
distortion  and  the  relaxed  structure  contains  two  four- 
membered  rings  per  unit  cell.  If  the  pairs  of  near¬ 
neighbor  atoms  are  assumed  to  be  bonded,  all  the  atoms 
have  fourfold  coordination  with  maximum  bond  length 
13%  greater  than  the  bulk  value.  Unlike  the  15  CSL 
state,  this  translation  state  does  allow  new  bonds  to  form 
on  both  sides  of  the  interface. 

Figures  1(e)  and  1(f)  illustrate  the  geometric  struc¬ 
tures  obtained  after  the  relaxing  of  two  15*  boundaries. 
Figure  1(e)  shows  the  CSL  after  relaxation.  VVe  find 
that  the  distortion  has  produced  three  atoms  per  unit  cell 
which  have  only  three  near  neighbors  and  the  structure 
contains  four  four-membered  rings  per  unit  cell.  The 
coincidence-site  atom  in  this  structure  has  five  near 
neighbors  as  a  result  of  formation  of  a  dimer  bond  with 
an  atom  in  the  same  layer  which  becomes  fourfold  coor¬ 
dinated.  The  occurrence  of  a  fivefold-coordinated  atom 
in  this  sfucture  is  particularly  interesting  in  the  context 
ol  recent  suggestions  that  such  defects  may  dominate 
over  threefold-coordinated  atoms  in  amorphous  silicon."' 
Finally,  Fig.  1(f)  shows  the  relaxed  15*  1^(210)  state. 
This  boundary  relaxes  by  forming  dimer  bonds  between 
atoms  in  the  layers  above  and  below  the  interface.  The 
dimers  bond  across  the  interface  to  form  continuous 
chains  of  fourfold-coordinated  atoms  lying  along  (310) 
as  can  be  deduced  from  the  figure.  However,  after 
reconstruction,  there  are  still  two  atoms  per  unit  cell  that 
have  only  three  near  neighbors  and  only  one  weak  four- 
membered  ring.  We  note  in  particular  from  Table  I  that 
the  boundary  undergoes  a  marked  contraction  in  width 
during  relaxation. 

For  both  the  15  and  15*  boundaries,  the  nonzero- 
translation  states  studied  had  the  lowest  total  energies. 
It  has  often  been  assumed  that  the  presence  of  the  coin¬ 
cidence  atom  m  a  site  belonging  to  both  crystals  would 
lead  to  a  low-energy  boundary.  Previous  experimental 
and  theoretical  results^  for  tilt  boundaries  in  semicon¬ 
ductors  have  shown  that  this  is  not  necessarily  the  case 
and  the  results  of  our  calculations  for  twist  boundaries  in 
germanium  support  this.  In  addition,  whereas  the  CSL 


states  both  involve  a  slicht  expansion  of  the  local  gram- 
boundary  volume,  the  translated  states  both  undergo  a 
contraction.  Previous  theoretical  studies  of  grain  bound¬ 
aries  in  metals  and  ionic  materials  have  always  predicted 
an  expansion  at  the  grain  boundary.  However,  small 
contractions  have  been  calculated  for  tilt  boundaries  in 
semiconductors  with  use  of  simple  valence  force-field  cal¬ 
culations"'  and  a  possible  small  contraction  in  a  22° 
(001)  twist  boundary  in  germanium  has  been  detected 
experimentally  by  use  of  electron  diffraction,  -'  The  mag¬ 
nitude  of  this  contraction  is  very  similar  to  that  reported 
here. 

In  any  simulated  quench  or  stimulated  anneal  per¬ 
formed  at  finite  cooling  rate  there  is  a  possibility  that  the 
system  will  not  reach  the  ground  state  but  will  instead 
remain  in  a  metastable  state  of  higher  energy.  To  test 
whether  our  lowest-energy  boundary  was  a  metastable 
state  we  made  a  large  perturbation  to  the  structure.  The 
perturbation  consisted  of  moving  all  the  atoms  in  the 
four  layers  around  the  interface  by  distances  of  10%  of 
the  bulk  bond  length.  The  perturbed  configuration  was 
then  allowed  to  relax  to  equilibrium.  Of  the  twenty  ions 
that  were  displaced,  eighteen  moved  closer  to  their  posi¬ 
tions  in  the  original  boundary:  but  two  atoms  in  the  layer 
immediately  below  the  interface  moved  further  from 
their  original  positions  in  a  direction  away  from  the  in¬ 
terface.  so  that  the  lengths  of  their  bonds  across  the  in¬ 
terface  were  increased  but  the  bonds  to  the  atoms  two 
layers  from  the  interface  were  reduced  in  length  The 
connectivities  of  the  two  relaxed  structures  were  identi¬ 
cal  and  the  boundary  formed  after  the  perturbation  was 
only  0.1  eV  higher  in  energy  than  the  original  boundary. 
This  then  leads  to  a  very  interesting  and  unexpected  re¬ 
sult.  The  structures  with  the  connectivity  of  the  original 
boundary  lie  m  an  "attractive  "  region  of  phase  space 
since  the  connectivity  is  not  destroyed  by  a  significant 
perturbation  to  the  system.  Moreover,  within  this  region 
of  phase  space  there  are  at  least  two  known  local  energy- 
minima  and  possibly  many  more.  Tunneling  transitions 
between  configurations  in  the  energy  minima  should  then 
give  a  range  of  effects  similar  to  those  produced  by  two- 
level  systems  in  glasses."'  Further  investigation  of  the 
energies  of  metastable  minima  and  tunneling  barriers  in 
this  region  of  phase  space  is  currently  underway.  Final¬ 
ly,  the  energy  difference  between  the  two  local  energy 
minima  is  smaller  than  the  energy  differences  between 
different  translation  states  and  this  suggests  that  the  15 
ST -'210)  translation  state  is  the  lowest-energy  boundary 
of  the  four  studied  here.  V\'ithout  calculating  the  ener¬ 
gies  of  many  other  translation  states,  it  is  not  certain 
that  the  ;ij'21(),'  translation  state  is  the  minimum- 
energy  boundary  However,  the  complete  fourfold  coor¬ 
dination  present  in  the  relaxed  structure  of  this  boundary 
must  make  it  a  good  candidate. 

Finally,  we  note  that  the  15  and  S5*  grain  boundaries 
are  related  to  one  another  hv  the  removal  or  insertion  of 
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an  odd  number  of  atomic  layers  normal  to  the  interface. 
To  fabricate  a  15  boundary  experimentally,  two  oriented 
crystals  with  (001)  surfaces  would  be  rotated  through 
36.9°  and  then  sintered  together  (see,  e.g..  Ref.  5).  In 
general,  the  resulting  boundary  will  be  a  mixture  of 
36.9°  and  —53.1°  rotations  unless  the  original  crystals 
had  monolayer-smooth  (001)  surfaces.  Consequently,  a 
fabricated  grain  boundary  would  normally  contain  both 
15  and  15*  interfaces. 
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